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ABSTRACT 

We report on millisecond variability from the X-ray transient XTE J1739-285. We detected six X-ray type 
I bursts and found evidence for oscillations at 1122 ± 0.3 Hz in the brightest X-ray burst. Taking into con- 
sideration the power in the oscillations and the number of trials in the search, the detection is significant at 
the 99.96% confidence level. If the oscillations are confirmed, the oscillation frequency would suggest that 
XTE J 1739-285 contains the fastest rotating neutron star yet found. We also found millisecond quasiperiodic 
oscillations in the persistent emission with frequencies ranging from 757 Hz to 862 Hz. Using the brightest 
burst, we derive an upper limit on the source distance of about 10.6 kpc. 

Subject headings: accretion, accretion disks — gravitation — relativity — stars: individual (XTE J1739-285) 
— stars: neutron — X-rays: stars 



1. INTRODUCTION 

Weakly magnetized neutron stars can be spun up to 
rates of seve ral 100 Hz by a ccretion in low-mass X- 
ray binaries dAlpar et ail 119821) . The first direct mea- 
surements of millisecond spin rates in actively accret- 
ing neutron stars in low-mass X-ray binaries (LMXBs) 
came from the discovery of oscillations in thermonu- 
clear X-ray bursts occ urring on the neutron star surface 
(IStrohmayer et al.l Il996h . The burst oscillation frequencies 
were later found to be nearly equal to those of coherent 
pulsa t ions din 't Zand et all |200U fStrohmaver & Markvvardtl 
I2002t IChak rabarty et al. 2003r iStrohmayer et al.l 120031) im- 
plying that the burst oscillation frequency indicates the neu- 
tron star spin rate. This X-ray technique has no known bi- 
ases against the detection of very high spin rates, unlike radio 
pulsation searches, and the sample of X-ray burst oscillation 
frequencies has been exploited to constrain the neutron star 
spin rate distribution. Analysis of a sample of 1 1 X-ray mea- 
sured spin frequencies in the range 270-619 Hz suggests a 
limiting spin rate near 760 Hz if the distribution is uniform 
and bounded (Chakrabartv et al. 2003). This is below the ex- 
pected maximum spin frequency possible without centrifugal 
breakup and has been interpreted as evidence that some phys- 
ical process, possibly gravitational radiation, limits the max- 
imum possible spin rate. However, the discovery of a radio 
pulsar spinning at 716 Hz, a frequency above any previously 
measured in X-rays, suggests that the true maximum spin rate 
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is higher dHessels et al.ll2006l) . 

Here, we describe observations made with the Rossi X- 
Ray Timing Explorer (RXTE; Bradt, Rothschild, & Swank 
1993) following the detection of X-ray bursts from the tran- 
sient source XTE J1739-285 as part of a program to search for 
millisecond oscillations in both X-ray bursts and persistent 
emission from neutron star X-ray binaries which are newly 
disco vered or found to be active (Ka aret et al.1 12002L 120031 
2006). We detected six X-ray bursts and found oscillations 
at a frequency of 1122 Hz in the brightest burst. This sug- 
gests that XTE J 1739-285 contains the most rapidly rotating 
neutron star yet discovered. We describe our observations in 
§2, the X-ray bursts in §3, and the persistent emission and the 
discovery of kHz QPOs from the source in §4. We discuss the 
results in §5. 

2. OBSERVATIONS OF XTE J1739-285 

XTE J 1739-285 is a transient neutron-star low-mass X- 
ray binary (NS-LMXB) that was discovered during RXTE 
PCA scans of the Ga lactic bulge on 1999 October 19 
dMarkwardt et al.l 1 19991) and underwent short outbursts in 
May 2001 and O ctober 2003. The sourc e became active again 
in August 2005 (Bodaghee et al. 2005) and two X-ray bursts 
were detected with the JEM-X inst rument on IN TEGRAL on 
2005 September 30 and October 4 dBrandt et al. 1120051) . 

Triggered by the detection of the X-ray bursts, we obtained 
19 observations using RXTE in the period beginning 2005 
October 12 and ending 2005 November 16. Data were ob- 
tained with the Proportional Counter Array (PCA) in a spec- 
tral mode (Standard 2) with 256 energy channels and 16 s 
time resolution, a low -resolution timing mode (Standard 1) 
with no energy information and 0. 125 s time resolution, and a 
high-resolution timing mode (Even t mode) with 122 u s time 
resolution and 64 energy channels dJahod a et al. 2006). 

3. X-RAY BURSTS 

We searched the Standard 1 data for X-ray bursts and found 
six, see Table [TJ We examined their evolution by extract- 
ing spectra for 0.25 s intervals of event mode data using all 
Proportional Counter Units (PCUs) that were on during each 
burst and all layers. The fluxes were corrected for dead time 
effects with a maximum correction of 5.5%. To remove the 



2 



Kaaret et al. 



TABLE 1 
Properties of X-Ray Bursts 



# Time Peak flux Decay 



1 Oct 31 07:59:09 1.6 ±0.3 5.0 

2 Nov 4 11:34:23 2.8 ± 0.5 5.4 

3 Nov 7 05:34:04 2.2 ± 0.4 5.0 

4 Nov 7 07:30:53 1.1 ± 0.3 3.0 

5 Nov 8 08:19:08 1.0 ±0.2 3.1 

6 Nov 11 09:47:10 1.2 ± 0.3 5.7 



NOTE. — The table lists for each burst: Time - the time (UTC) at the start of each 
burst (all bursts occurred in 2005); Peak flux - the bolometric peak flux corrected for 
absorption in units of 10 ~ 8 erg cm -2 s~ 1 ; and Decay - the decay (e-folding) time in 
seconds. 

contribution of the persistent emission, we subtracted off a 
spectrum from 10 s of data preceding each burst. Spectra 
were fitted in the 3-18 keV band (channels 3-24 in the event 
mode data) with an absorbed blackbody model with the col- 
umn density fixed to Nh = 7.5 x 10 21 cm -2 which is the 
hydro gen column density along the line of sight in the Milky 
Way dDickev & Lock man 1990). We found a burst-like event 
at Nov 11 10:27:22 UTC which appears in only one of four 
PCUs on at the time and is likely a detector breakdown event 
and not an X-ray burst. 

The brightest burst is # 2, which had a peak bolomet- 
ric flux, as measured in 0.25 s intervals, of (2.8 ± 0.5) x 
10 -8 erg cm~ 2 s _1 . The blackbody temperature rose rapidly 
during the flux rise, reached a maximum near 2.5 keV, and 
then decayed as the flux decayed. This behavior is indicative 
of heating during the rise and cooling during the decay and 
is characteristic of type I X-ray bursts. iKuulkers et alj (2003) 
found that the empirical maximum bolometric peak luminosi- 
ties for photospheric radius expansion bursts from a sample of 
12 bursters located in globular clusters with known distances 
was 3.8 x 10 38 ergs -1 , with an accuracy of ~15%. None of 
our bursts showed evidence of photospheric radius expansion 
and, thus, their luminosities should be below this value. The 
peak flux of the brightest burst implies an upper limit on the 
distance to XTE J1739-285 of about 10.6 kpc. This is consis- 
tent w ith, but somewhat closer than, the upper limit of 12 kpc 
from lTorres et all (120061) . 

Bursts 1 and 4-6 have similar peak fluxes, while bursts 2 
and 3 were brighter, see Table Q] Burst 3 had a rise time near 
3 s, while the other rise times were near 1 s. Burst 6 had a 
relatively broad maximum extending over 5 s. The decay (e- 
folding) times are in the range 3-6 s. The fact that the bursts 
have fast rises and durations of tens of seconds suggests that 
they are typical helium bursts. 

Using high time resolution data, we computed power spec- 
tra for overlapping 4 s intervals with 0.125 s between the 
starts of successive intervals using events in PHA channels 
4-24 for PCU0 and channels 3-24 all other PCUs. The back- 
ground, particularly at low energies, is higher in PCU0 due 
to the loss of its propane layer. We searched over the du- 
ration of each burst defined as the interval where the total 
counting rate exceeds the background plus persistent count- 
ing rate, measured a few seconds before the burst, by at least 
a factor of 2. We searched for excess power in the range 20- 
2000 Hz, corresponding to 7920 frequency bins in the 4 s 
FFTs. We found osc illations in the s econd burst with a Leahy 
normalized power (Leahy et al. 1983) of 42.82 at a frequency 
of 1122 ± 0.3 Hz, occurring in the burst decay 15 s after 
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FIG. 1 . — Dynamical power spectrum of burst 2 with the burst light curve 
(count rate) and the power at 1122.0 Hz superposed. The contours are at 
Leahy powers of 12, 20, 28, and 36. The contours are generated from power 
spectra for overlapping 4 s intervals of data with the power plotted at the 
midpoint of the interval. The curve which peaks near 2 s is the rate in the 3- 
18 keV band from all the PCUs on during the observation (PCU0 and PCU2). 
The curve which peaks near 15 s is the Leahy power at 1122.0 Hz. The 
scales for both curves are on the right side. The oscillation is significant at 
the 99.96% confidence level. 

the burst rise, see Fig. Q] The four consecutive FFTs around 
the one with the peak power all have powers above 40.31 at 
1122.0 Hz. 

To evaluate the significance of the signal, we used a Monte 
Carlo simulation which generates Poisson distributed events 
following the light curve from burst 2 in 0.125 s bins after 
smoothing with a moving average over 9 bins. The deadtime 
of the PC A is modeled by removin g any event which occurs 
within 10 fis after a previous event dlahoda et ai]|2006f) . The 
number of events generated in each time bin is larger than 
the observed counts so that after the deadtime correction the 
number matches that in the actual light curve within Poisson 
fluctuations. We generated 400,000 trial bursts. Each simu- 
lated burst was analyzed using exactly the same analysis done 
on the real data, specifically by calculating 4 s FFTs at over- 
lapping intervals with starts each 0.125 seconds and searching 
for power in the 20-2000 Hz interval. The chance probability 
of occurrence of the observed signal is calculated by counting 
the fraction of trial bursts with powers equal to or exceed- 
ing those observed. This procedure eliminates the ambiguity 
in estimating the equivalent number of independent trials for 
the overlapping FFTs. We found 58 bursts with at least one 
power above 42.82 of which 14 had at least 4 consecutive 
powers above 40.31 at the same frequency. We estimate the 
chance probability of occurrence of the observed signal to be 
3.5 x 10~ 5 . Allowing two trials for two energy bands, the 
probability is 7.0 x 10~ 5 equivalent to a 3.97er confidence 
level. 

The signal at 1122 Hz was detected in the brightest burst. 
Since the bursts are not a uniform sample and the presence 
of oscillations may depend on burst properties and/or accre- 
tion history (e.g. Watts, Strohmayer, Markwardt 2005), one 
can reasonably argue that the significance should be evalu- 
ated for each burst individually, particularly in this case since 
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TABLE 2 

High frequency peaks in the persistent emission 



10 15 20 
Time [s] 

FIG. 2. — Power at 1122 Hz in independent 1 s intervals (crosses). The 
solid line is the count rate in the same 1 s intervals. The powers at 1122 Hz 
used in the significance calculation are, in time order starting at 12.5 s, 3.90, 
9.51, 9.34, 15.37, and 9.67. We note that the power at 17.5 s at 1123 Hz is 
not shown on the plot and is 7.94. 

the brightest burst is the most likely to give a detectable sig- 
nal. However, a more conservative approach is to consider 
the full set of six bursts. The durations of the bursts, as de- 
fined above, are 21, 21, 22, 10, 11, and 18 s for bursts 1-6, 
respectively. To account for all the trials in all the bursts, we 
multiply the number of trials by an additional factor of 103/21 
= 4.9 for a chance probability of 3.4 x 10~ 4 . The oscillation 
is detected at a confidence level of 99.966%, equivalent to a 
3.6a significance. 

To obtain an estimate of the significance of the oscillations 
which does not depend on a simulation, we calculated FFTs 
in independent 4 s, 2 s, and 1 s intervals of data from burst 2. 
There are five successive 1 s FFTs with power at 1 122 Hz, see 
Fig. |2 The single trial probability for having five powers at 
or above the observed levels in the 1 s FFTs is 4.2 x 10~ n . 
Accounting for all sequential combinations of FFTs of fixed 
length of 1, 2, or 4 s with a total duration of 5 s or less cov- 
ering the 20-2000 Hz range in two energy bands, we esti- 
mate a total of 605880 trials: (17+18+19+20+21)x 1980x2 + 
(9+ 10)x 3960x2 + 5x7920x2. The chance probability tak- 
ing into account the number of trials in the one burst is then 
2.5 x 10 -5 (4.2cr). Accounting for the other bursts, the chance 
probability is 1.2 x 10~ 4 (3.8c). This probability agrees 
within a factor of 3 with that in the previous paragraph. 

In the 4 s interval with the maximum power, the fractional 
rms amplitude of the 1122.0 Hz oscillation is 0.13 ± 0.02. 
The quality factor, Q = ty/Av, of the oscillation in burst 2 
is Q > 1000. This is consistent with the values found in 
burst oscillations and inconsistent with the values found in 
high frequency QPOs in the persistent emission. We searched 
for oscillations near half the frequency. The strongest signal 
has a Leahy power of 13.6 at 561.25 Hz and occurs 7.75 s 
after the burst rise. This is not a significant signal; the chance 
probability of occurrence within the narrow window searched, 
562.5-558.5 Hz, is 0.09. At the time of the maximum strength 
of the 1122 Hz signal, the Leahy power near 561 Hz is less 
than 6. 



Time 


Centroid 


Width 


Amplitude 


(UTC) 


(Hz) 


(Hz) 


(%) 


Oct 31 06:22:11 


837.9±0.3 


5.3±0.9 


1.7±0.2 


Nov 01 02:48:14 


784.6±1.1 


14±3 


2.4±0.4 


Nov 08 04:37:47 


814.4±3.2 


36±9 


3.1±0.6 


Nov 08 09:46:31 


756.7±0.7 


9±3 


2.1 ±0.4 


Nov 10 08:30:35 


861.8±0.7 


13±2 


2.7±0.3 


Nov 10 10:04:43 


846.0±2.9 


33±8 


2.7±0.5 



NOTE. — The table includes: Time - the UTC time at the beginning of the obser- 
vation, all observations were in 2005; Centroid and Width - of the fitted Lorentzian; 
Amplitude - RMS fractional amplitude of the fitted Lorentzian. 



4. PERSISTENT EMISSION 

To study the persistent emission spectra, we used Standard- 
2 data from PCU 2, removing data around the X-ray bursts, 
and estimating the background using bright source back- 
ground files. The source was in a relatively hard state, with 
a hard X-ray color, HC, defined as the count rate in the 9.7- 
16 keV band divided by the 6.0-9.7 keV rate, of HC > 0.35 
for all observations before October 25. The flux was roughly 
constant and near 4 x 10~ 10 erg cm" 2 s" 1 in the 2-20 keV 
band. The spectra were adequately described by the absorbed 
(Nh = 7.5 x 10 21 cm -2 ) sum of a power-law with a pho- 
ton index near 2 and a gaussian emission line with a centroid 
of 6.7 keV and a flux consisten t with that from Galactic ridge 
emiss ion at the source position ( Rev nivtsev. Molkov. Sazonovl 
2006). After October 25, the source stayed in a softer state, 
HC < 0.35. The flux varied by a factor of two, with the maxi- 
mum near 1.4 x 10~ 9 erg cm -2 s" 1 in the 2-20 keV band. The 
spectra showed distinct curvature at high energies and were 
adequately fitted with the absorbed sum of a Comptonization 
model (compst in xspec) with a temperature in the range 4— 
12 keV and an optical depth in the range 3-10, a blackbody 
with a temperature near 1 .6 ke V, and a gaussian emission line 
with a centroid of 6.7 keV, a width less than 0.9 keV, and 
an equivalent width of 100-300 eV, in some cases larger than 
the expected Galactic ridge emission. The blackbody compo- 
nent was not required for some lower flux observations. The 
spectra a re similar to those of othe r moderate luminosity NS- 
LMXBs dKaaret et alJ200ll200l . 

To further investigate the source state, we produced low fre- 
quency power spectra using event mode data from all PCUs 
on during each observation. The sum power spectrum for ob- 
servations with HC > 0.35 shows very-low-frequency noise 
(VLFN) and so-called high-frequency noise (HFN). The to- 
tal noise power in the 0.01-100 Hz band is 1.0% fractional 
rms. We fit the power spectrum with a model consisting of a 
power-law to describe th e VLFN and an exp onentially cutoff 
power-law for the HFN (Ivan der Klislll995l) . The best fit pa- 
rameters are a power law index for the VLFN of 1.81 ± 0.13 
and a power law index of 0.8 ± 0.1 and a cutoff frequency in 
the range of 20-60 Hz for the HFN. The weakness of the sig- 
nal at high frequencies prevents an accurate determination of 
the HFN parameters, but a HFN noise component is required, 
A\ 2 — 49. This power spectrum is consistent with those in 
the "lower banana" state. There is no detectable HFN in the 
power spectrum from the observations with the softer hard 
color, HC < 0.35. The total noise power in the 0.01-100 Hz 
band is 0.2% fractional rms and the spectrum is adequately 
described by a single power law with an index of 1.08 ± 0.02. 
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FIG. 3. — Power spectrum showing a kHz QPO from October 31. 
power is Leahy normalized. 
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The shape of the power spectrum is consistent with those seen 
in the "upper banana" state, but the total noise power is lower 
than usual and s pectrally the state would be classified as the 
"lower banana" (Ivan der Klisl[l995l) . 

Based on the timing and color information, we identify 
XTE J1739-285 as an atoll source. The source appears to have 
been in the "banana" state during all of the observations ana- 
lyzed here, although its evolution along the "banana" appears 
somewhat unusual. 

We searched for high frequency quasiperiodic oscillations 
(QPOs) in each uninterrupted RXTE observation window. We 
calculated averages of 2 s power spectra for PCA events in the 
4.2-20.4 keV energy band (channels 6-26 in the event mode 
data). We searched for peaks in the range from 100 Hz to 
3000 Hz and fit any peak found with a Lorentzian plus a con- 
stant equal to the calculated the Poisson noise level. There 
were several strong QPO detections, see Table [2] The sec- 
ond most significant detection occurred on October 3 1 and is 
shown in Fig. [3] Allowing for 1 102 trials, which is calculated 
by dividing the search interval of 2900 Hz by the QPO trial 
widths of 5, 10, 20, 50, and 100 Hz, we estimate a chance 
probability of occurrence of 1.2 x 10~ 13 for this QPO. The 
chance probabilities of occurrence, taking into account the 



number of trials, of all the other QPOs listed in the Table are 
less than 1 x 10~ 4 . These QPO detections establish XTE 
J 1739-285 as a new member of the class of neutron-star low- 
mass X-ray binaries producing kHz QPOs. We did not detect 
two simultaneous kHz QPOs in any observation. 

5. DISCUSSION 

These observations establish that XTE J1739-285 exhibits 
millisecond oscillations in its persistent X-ray emission. The 
properties of the source are generally similar to those of other 
atoll neutron-star X-ray binaries, although the evolution of the 
timing noise in the banana state is unusual. 

The signal at 1122 Hz present in burst 2 is significant at 
the 99.96% confidence level. If this signal represents a true 
burst oscillation then it would be of substantial interest. The 
near equality of the burst oscillation frequency with the fre- 
quency of coh erent pulsations in the millisecond pulsars SAX 
J1808.4-3658 din 't Zand et al.l200Ulc¥afaabartv et alj2003l) 
and XTE J1814-338 dStrohmaver et alj 120031) and the fre- 
quen cy of coherent pulsations in a s uperburst from 4U 1636- 
536 dStrohmaver & M arkwardt 2002) strongly suggests that 
the burst oscillation frequency indicates the neutron star spin 
frequency. The lack of any significant signal near 561 Hz in 
the burst from XTE J 1739-285 supports this interpretation and 
suggests that the possible 1 122 Hz oscillation would be most 
naturally interpreted as the spin rate of the neutron star. 

If the burst oscillation frequency of 1 122 Hz is the spin rate 
of the neutron star, then XTE J1739-285 contains the most 
rapidly rotating neutron star yet discovered. This spin rate 
is close to the centrifugal breakup limit for some equations 
of state of nuclear matter (Burgi oTSchulze. & Weberl 120031) 
and, therefore, may remove the motivation for a physical limit 
on neutron star spin other than the centrifugal breakup limit. 
Furthermore, such a high spin rate would place constraints on 
the nuclear equation of state, particularly if combined with a 
measurement of the mass and/or radius of the neutron star. 
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